Jahn- Teller like origin of the tetragonal distortion 
in disordered Fe-Pd magnetic shape memory alloys 
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The electronic structure and magnetic properties of disordered Fej^Pdioo-i^ alloys (50 < a; < 
85) are investigated in the framework of density functional theory using the full potential local 
orbital method (FPLO). Disorder is treated in the coherent potential approximation (CPA). Our 
calculations explain the experimental magnetization data. The origin of the tetragonal distortion 
in the Fe-Pd magnetic shape memory alloys is found to be a Jahn- Teller like effect which allows the 
system to reduce its band energy in a narrow composition range. Prospects for an optimization of 
the alloys' properties by adding third elements are discussed. 



Magnetic shape memory (MSM) alloys have attracted 
considerable attention as materials for actuator and 
sensor applications, due to large magnetically induced 
strains of up to 9.5% in Ni2MnGaJ. A promising MSM 
alloy is disordered FeroPdso^ with a relatively high block- 
ing stress, a high saturation magnetization and a Curie 
temperature well above room temperature, as well as a 
high ductility, in contrast to Ni2MnGa. To date, mag- 
netic field induced strains of about 3% have been ob- 
served in Fe-Pd alloys^ 

The unusually large strains observed in Ni2MnGa and 
Fe7oPd3o, which are one to two orders of magnitude 
higher than in conventional magnetostrictive materials, 
arc attributed to a reorientation of martensitic twins in 
the magnetic field. Twin variants which have their easy 
axis of magnetization aligned along the field direction 
gain Zeeman energy and grow at the cost of variants 
with an unfavorable alignment. This effect thus requires 
a magnetocrystalline anisotropy energy large enough to 
prevent a rotation of the magnetization within a twin 
variant and a high mobility of twin boundaries. 

In this Letter, we investigate the electronic structure of 
disordered Fea;Pdioo-K alloys (50 < a; < 85) by means of 
density functional theory calculations. In particular, we 
explore the electronic origin of the tetragonal distortion 
resulting in the fct martensite phase, where the MSM 
effect is observed. On the basis of these results, we dis- 
cuss prospects for a further optimization of the alloys' 
properties. 

The value of the c/a-ratio in the fct martensite phase 
determines the MSM properties of Fe-Pd alloys in a 
twofold way: First, the maximum strain is given by 
e = Al/l = |1 — c/a\ under the assumption of a hun- 
dred percent conversion of twin variants. Second, the 
magnetocrystalline anisotropy energy, which determines 
the blocking stress, is expected to depend strongly on the 
c/a-ratio4i^ Thus, understanding the mechanism behind 
the tetragonal distortion is important for the develop- 
ment of alloys with improved MSM properties. 

The martensitic transformation from fee to fct in dis- 
ordered Fe-Pd alloys is observed in a narrow composition 
range between 29 and 32 at% Pdr^ close to the transi- 



tion from an fee to a bcc ground state with increasing Fe 
content. Upon further cooling an irreversible transition 
from fct to a (bcc-like) bet phase occurs. The martensite 
temperature Tm is close to room temperature for opti- 
mal alloy composition, but depends strongly on the Fe 
content 1^ For applications a martensite temperature Tm 
clearly above room temperature would be desirable, since 
the operation range of MSM elements is limited by Tm- 
Attempts to raise Tm in the Fe-Pd system by addition of 
Co, Ni, Rh or Pt have so far been unsuccessfuljiiSiiSiiiS Fur- 
thermore, the c/a-ratio in the martensite phase depends 
both on composition and on temperature. A gradual in- 
crease of c/a from about 0.92 at -80° C to 0.97 near Tm 
was observedfi similar to the temperature dependence of 
c/a in FcsPt, but distinct from the prototype MSM alloy 
Ni2MnGa4 

We have calculated the electronic structure of disor- 
dered Fe^^Pdioo-z: alloys (50 < x < 85) in the frame- 
work of density functional theory (DFT) in the local 
spin density approximation (LSDA) using the high pre- 
cision full potential local orbital (FPLO, version 5.00-19) 
methodii in the scalar relativistic approximation. Disor- 
der was treated within the coherent potential approxima- 
tion (CPA)i^ and the LSDA parameterization of Perdew 
and WangH was employed. The basis set consisted of Fe 
3spd4sp and Pd 4spd5sp optimized local orbitals. The 
exponent of the confining potential in the basis optimiza- 
tion was set to n = 6, which provided lower total energies 
and a better numerical stability compared to the default 
value n = 4. For the k-space integrations 1728 points 
in the full Brillouin zone (FBZ) were used, and the re- 
sults were cross-checked with 8000 points in the FBZ. 
Additionally, we have performed calculations with the 
Akai-KKR codei^ in the muffin tin approximation using 
the generalized gradient approximation (GGA) to DFT 
in the parameterization by Perdew et alJ^ 

The calculated magnetic moments of fee Fe^^Pdioo-^: 
alloys (50 < x < 70) are in an overall good agreement 
with the experimental data by Matsui et ali^ The aver- 
age spin moment per atom, calculated with FPLO at the 
respective experimental lattice parameter, increases ap- 
proximately linear with increasing Fe content from 1.55 
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FIG. 1: Total energy as a function of c/a-ratio for Fe8oPd2o 
calculated with the atomic volume fixed to Vat = 13.256 A"^. 
The total energy exhibits a local minimum at c/o=0.875 cor- 
responding to an fct martensite structure. 



Hb for a; = 50 to 1.98 foi' x ~ 70. Adding the con- 
tributions of the orbital moments, which are estimated 
from relativistic supercell calculations^ to be about 0.13 
jjLB for Fe and 0.04 jiB for Pd, the calculated moments 
agree within 0.05 hb with the experimental values for 
the average total moments 1.61 ^b {x = 50) and 2.04 ^b 
{x = 70). 

Assuming that the individual atomic moments do not 
depend on the compostion, Matsui et alt^ had to use 
an unusually large Pd moment of 0.56 /is to fit the con- 
centration dependence of their magnetization data. This 
value is nearly twice the value wc obtained for the Pd spin 
moment (0.29 /is). Our calculations provide a different 
explanation for the observed concentration dependence. 
The spin moment of Fe decreases with increasing Fe con- 
tent due to magneto volume effects, from about 2.82 ^ib 
for X = 50 to 2.71 ^b for x = 70. A yet more pronounced 
reduction of the Fe moments with increasing Fe content 
was recently observed in neutron scattering experiments 
on Llo-ordered FePt alloys^i. 

The fee structure can be continuously transformed into 
a bcc structure by a change of the c/a-ratio, the so called 
Bain path. A (meta) stable fct phase must exhibit a 
(local) minimum in the free energy along the Bain path 
for 0.71 < c/a < 1. In our total energy calculations we 
find such a local minimum only in a narrow region of 
phase space, where the total energy differences between 
fee and bcc are small and the energy landscape along 
the Bain path is accordingly flat (Fig. [1]). This is in 
agreement with the experimentally observed proximity of 
the martensite phase to the transition from an fee to a bcc 
ground state. While in the LSDA calculations the fcc-bcc 
transition is found at a somewhat higher Fe concentration 
X « 80 in the vicinity of the experimental volume, the 
KKR GGA data yield almost exactly the experimental 
value x ~ 70. 

Further insight into the origin of the martensite phase 
is obtained from the density of states (DOS) shown in 
Fig. O The DOS in the vicinity of the Fermi energy is 
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FIG. 2: (Color online) DOS of Fe8oPd20 in the fee (sofid lines) 
and fct (dashed lines) structures corresponding to Fig. [T] The 
total DOS (thick lines) and the partial Fe 3d DOS (thin lines) 
are shown. The minority spin DOS (J.) exhibits a peak at the 
Fermi energy (vertical line) in fee which splits in fct. 



dominated by Fe 3d states, which arc hybridized with 
Pd 4d states (not shown). The majority spin channel 
of the Fe 3d subshcU is completely filled, while in the 
partially filled minority spin channel a peak is located 
right at the Fermi energy in the fee austenite phase. In 
the fct martensite phase this peak is split and weight 
is shifted below and above the Fermi energy, indicative 
of a cooperative Jahn- Teller effect. The gain in band 
energy due to this splitting (estimated from a single step 
calculation) is a few tenths of a mHartree. 

Further support for a Jahn- Teller like origin of the 
tetragonal distortion in the disordered Fe-Pd MSM al- 
loys is obtained from the spectral density of the minority 
spin states shown in Fig. [3] For the cubic austenite phase 
(upper panel) the threefold degenerate t2g states of the Fe 
3d minority spin shell are located right at the Fermi level 
close to the F-point (note that due to disorder a broad- 
ening of the bands occurs). The tetragonal distortion in 
the martensite phase (lower panel) causes a splitting of 
the t2g states, which are shifted away from the Fermi en- 
ergy, resulting in a gain in band energy in a Jahn- Teller 
like fashion. 

On the basis of the results discussed above, the follow- 
ing conclusions can be drawn: The presence of degen- 
erate Fe 3d minority spin bands at the Fermi energy in 
Fe-Pd MSM alloys favors a structural distortion, which 
is stabilized by a gain in band energy. However, due to 
band dispersion and disorder, this energy gain is small 
and thus can only compete with other energy terms (like 
the Madelung energy favoring close packed structures) if 
the energy landscape is sufficiently fiat, i.e. close to the 
fcc-bcc transition. The gain in band energy due to the 
Jahn- Teller like splitting is approximately proportional 
to the height of the peak in the minority spin DOS. The 
energy difference between fee and bcc changes by approx- 
imately 4 meV/at%, while the height of the peak in the 
DOS shows only a moderate concentration dependence. 
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FIG. 3: CPA minority spin spectral densities of FegoPd2o in 
the fee structure (upper panel) and fct (lower panel) corre- 
sponding to Fig. [1] 



Thus already small deviations from the optimal stoi- 
chometry result in relatively steep energy curves along 
the Bain path and cause the local fct minimum to disap- 
pear. Consequently, a stabilization of the fct structure 
in terms of c/a-ratio and martensite temperature Tm is 
expected with increasing iron concentration x, up to the 
point where a bee (or bet) structure becomes favorable. 
Alternatively, alloying Fe-Pd with ternary elements could 
lead to a stabilization of the fct phase. According to our 
KKR GGA calculations, Fe-Pd-X alloys with about 5 
at% of the ternary elements X= W, Au or Pt show a 
more pronounced peak in the minority spin DOS close to 
the respective fcc-bcc transition and are therefore inter- 
esting candidates for MSM alloys with enhanced T^- 

Our results imply that Tm and c/a in Fe-Pd(-X) alloys 
are determined by a delicate balance between the energy 
gain due to the Jahn- Teller splitting and the proximity to 
the fcc-bcc transition. This is in line with the experimen- 
tally observed T„i of binary and ternary Fe-Pd(-X) MSM 
alloys (X= Co, Ni, Pt), which shows a strong increase 
with the Fe concentration x for alloys having a compa- 
rable amount of the ternary element X (Fig. [4]). The 
critical Fe concentration Xc, where fee and bcc are de- 



generate in energy depends on the amount and the type 
of X (Fig. 31). The values of Xc have been obtained from 
KKR GGA calculations and depend below 10 at% of X 
approximately linear on the amount of ternary element. 
The KKR GGA calculations yield Xc for the binary al- 
loys in excellent agreement with the experimental trends, 
and the relative trends upon alloying can be considered 
trustworthy despite the limited accuracy of the muffin tin 
approximation. For Co and Ni, Xc is shifted to lower val- 
ues, while for Pt it is increased with respect to the binary 
system. This naturally explains that in the experiments 
of Ref. [13 addition of Pt at constant Fe content was found 
to have a negative influence on T^, as with increasing Pt 
content Xc was shifted away from 70 at% Fe concentra- 
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FIG. 4: Experimental martensite temperature as a func- 
tion of the Fe concentration x for Fe^iPdioo-a; (a: Ref. [g), 
Fej;Pdioo-x-!/Nij, (b: Ref. H), Fe^Pdioo-x-i/Coy (c: Ref. H; 
d: Ref.i), and Fe^Pdioo-^-aPty (e: Ref.[l3, f: Ref.i). Num- 
bers near data points indicate the amount of ternary element 
y. Also shown is the calculated Fe concentration Xc for the 
fcc-bcc transition (dashed vertical lines) for Fe-Pd and ternary 
alloys with y= 5 at%. 



tion used. The somewhat higher value of Tm observed in 
Ref.i is in line with our argument. Consequently, a care- 
ful optimization of the Fe content is necessary to achieve 
Tm enhancement in ternary Fe-Pd-X alloys. 
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